Multi-User Cooperative Diversity through Network 
Coding Based on Classical Coding Theory 

Joao Luiz Rebelatto, Student Member, IEEE, Bartolomeu F. Uchoa-Filho, Member, IEEE, Yonghui Li, Senior 

Member, IEEE, and Branka Vucetic, Fellow, IEEE 



o 
in 



CO 



> 

in 
(N 

o 
o 



X 



Abstract — In this work, we propose and analyze a generalized 
construction of distributed network codes for a network consist- 
ing of M users sending different information to a common base 
station through independent block fading channels. The aim is 
to increase the diversity order of the system without reducing its 
throughput. The proposed scheme, called generalized dynamic- 
network codes (GDNC), is a generalization of the dynamic- 
network codes (DNC) recently proposed by Xiao and Skoglund. 
The design of the network codes that maximize the diversity order 
is recognized as equivalent to the design of linear block codes over 
a nonbinary finite field under the Hamming metric. We prove that 
adopting a systematic generator matrix of a maximum distance 
separable block code over a sufficiently large finite field as the 
network transfer matrix is a sufficient condition for full diversity 
order under link failure model. The proposed generalization 
offers a much better tradeoff between rate and diversity order 
compared to the DNC. An outage probability analysis showing the 
improved performance is carried out, and computer simulations 
results are shown to agree with the analytical results. 

Index Terms — Cooperative communication, diversity order, 
nonbinary network coding, Singleton bound. 



I. Introduction 

IN a cooperative wireless communications system with 
multiple users transmitting independent information to a 
common base station (BS), besides broadcasting their own 
information, users help each other relaying their partners in- 
formation |T)-||7]. In the decode-and-forward (DAF) relaying 
protocol 0, the codeword relayed to the base station (BS) 
is a re-encoded version of the previously decoded codeword 
received in the broadcasting phase. A simple 2-user DAF 
network is illustrated in Fig.Q] In the first time slot (Fig. |l(a)) , 
each user broadcasts its own information (usually through 
orthogonal channels). In the second time slot (Fig. |l(b)| >, 
each user transmits its partner's information to the BS after 
decoding and re-encoding it. If the user cannot decode the 
partner's information, it sends its own information again. 
By transmitting the same information through independent 
channels, a cooperative diversity is achieved 0~|, 0. 

Network coding (9), iflOl , a method originally proposed to 
attain maximum information flow in a network, has recently 
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Fig. 1. Two-user decode-and-forward cooperative network, (a) Each user 
broadcasts its own information and (b) each user transmits its partner's 
information after decoding and re-encoding it. 



been applied to cooperative wireless communications systems 
to improve the bit error rate (BER) performance H-0, ifTTl - 
lfl4l . In a network coded system, relays process information 
from different users and perform linear combinations of the 
received signals, with coefficients chosen from a finite field 
GF( 9 ). 

A two-user cooperative system that employs binary network 
coding (BNC) was proposed in 0. In that scheme, each user 
transmits the binary sum (XOR) of its own source message and 
the message received from its partner (if correctly decoded). 
However, the scheme proposed in does not improve the 
system diversity order. 

In lfTTll . it was shown that binary network coding is not 
optimal for achieving full diversity in multiple user-multiple 
relay systems when link failures are taken into account. A sim- 
ilar result was shown in [5], however, instead of considering 
dedicated relays, therein the M users themselves act as relays 
for each other. The scheme proposed in 0, called dynamic- 
network codes (DNC), considers a fixed nonbinary network 
code. In the first time slot, each user broadcasts a single packet 
of its own to the BS as well as to the other users, which try 
to decode the packet. From the second time slot until the M th 
time slot, each user transmits to the BS M — 1 nonbinary linear 
combinations of the packets that it could successfully decode. 
With DNC, by using an appropriately designed network code, 
the diversity order was shown to be higher than that in binary 
network coded systems. The scheme is called "dynamic" in 
the sense that the network code is designed to perform well 
under a possible occurrence of errors in the channels between 
the users (inter-user channels). However, this design approach 
can become extremely complex as the number of users M 
grows, as shown in 0. 

In fl5l . lfl6l . Yeung and Cai derived a generalization of the 
Hamming, the Singleton and the Gilbert- Varshamov bounds, 
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extending concepts from the classical point-to-point coding 
theory to a communication network, and focusing on its 
error correction capability. In ifTTll , it was shown that the 
minimum distance of a network error-correction code for a 
single-source multicast case plays the same role as it plays in 
classical coding theory. The decoding performance under the 
occurrence of different types of errors and erasures was also 
evaluated. 

In this paper, we investigate another relationship between 
network codes and classical error-correcting codes. We elab- 
orate on the DNC scheme by first recognizing the problem 
as equivalent to that of designing linear block codes over 
GF(q) for erasure correction. In particular, for perfect inter- 
user channels, we note that the diversity order equals the 
minimum Hamming distance of the block code, so the network 
transfer matrix should correspond to the generator matrix of an 
optimal block code under the Hamming metric. The Singleton 
bound appears as a natural upper bound on the diversity 
order, and this bound is achieved with a sufficiently large field 
size Ifl8l . 

We then extend the DNC scheme by allowing each user to 
broadcast several (as opposed to just one) packets of its own in 
the broadcast phase, as well as to transmit several nonbinary 
linear combinations (of all correctly decoded packets) in the 
cooperative phase. From the block coding perspective, the so- 
called generalized dynamic-network codes (GDNC) consider a 
longer codeword, with more parity symbols, which improves 
the Singleton bound. 

Codes that achieve the Singleton bound are called maximum 
distance separable (MDS) codes [18]. Regarding the GDNC 
network code design, we show that if a generator matrix of a 
MDS code is used as the GDNC network code, the maximum 
diversity order is guaranteed. We also show that a much better 
tradeoff between rate and diversity order can be achieved, e.g., 
it is possible to achieve both rate and diversity order higher 
than those in the DNC scheme. 

The rest of this paper is organized as follows. The next 
section presents the system model and some relevant pre- 
vious works, including binary network coded cooperative 
systems []4j and the DNC scheme 0. The motivation for 
GDNC is presented in Section [HI] Section [TV] presents the 
proposed GDNC scheme, showing firstly a simple 2-user 
network and then generalizing it to M users. An outage 
probability analysis and our network code design approach are 
carried out in Sections [V] and [VTJ respectively. Simulations 
results are presented in Section IVIII Finally, Section IVIIII 
presents our conclusions and final comments. 

II. Preliminaries 

A. System Model 

The system consists of a multiple access channel (MAC) 
network, in which multiple users (M > 2) have different 
information to send to a common base station (BS). One 
time slot (TS) is defined as the time period in which all the 
M users perform a single transmission (through orthogonal 
channels, either in time, frequency or code). That is, one TS 
corresponds to M transmissions from individual users. The 



received baseband codeword at user i at time t is given by 

Vj,i,t = hjjjXjj + nj t ^t, (1) 

where j 6 {1, • • • , M} represents the transmit user index and 
i G {0, 1, • • • , M} the receive user index (0 corresponds to 
the BS). The index t denotes the time slot. Xj :t and yj,i,t are 
the transmitted and the received codewords, respectively, rij i t 
is the zero-mean additive white Gaussian noise with variance 
Nq/2 per dimension. The channel gain due to multipath fading 
is denoted by /ij,^, and it is assumed to have independent 
identically distributed (i.i.d.) (across space and time) Rayleigh 
distribution with unit variance. 

Assuming Xjj's to be i.i.d. Gaussian random variables and 
considering all the channels with the same average signal-to- 
noise ratio (SNR), the mutual information Ij t it between Xj t t 
and yj^t is 

/,. ,/ log(l • /',.,.» "SNR:. (2) 

In this work, we do not consider a joint network-channel 
code design Q, iTPH . [fT~3l - It is assumed that powerful enough 
channel codes are being used at a link level, so that xjj can 
be correctly decoded at user i if Ij,%,t > Pj.j.i, where rj ij is 
the information rate from User j to User i in the time slot t. 
Considering that all the users have the same rate, the index of 
r can be dropped. Thus, Xjj cannot be correctly decoded at 
user i if 

l%»,t| 2 <5, (3) 

where g = 2 SNR X . The probability that such an event happens is 
called the outage probability. For Rayleigh fading, the outage 
probability is calculated as 0, lfl9l 

P e = Pr { |%. M | 2 < <?} = 1 - e~ 9 « g, (4) 

where the approximation holds for a high SNR region. Con- 
sidering block fading, the diversity order D is defined as |fl9l 



SNR-+00 log SNR ' 

where P Q is the system overall outage probability. 

In this work, block fading means that fading coefficients are 
i.i.d. random variables for different blocks but constant during 
the same block. It is also assumed that the receivers have 
perfect channel state information (CSI), but the transmitters 
do not have any CSI. 

Before we continue, we establish that throughout this paper, 
the operators + and — represent operations over real numbers, 
© is the binary addition (XOR), and EE and EE are the addition 
and subtraction operations over a nonbinary field, respectively. 

B. Binary Network Coded Cooperation 

Instead of transmitting only the partner's information in 
the second TS, as shown in Fig. Q] each user of this 2-user 
system can transmit a binary sum of its own information and 
its partner's information. Considering for the moment that the 
two inter-user channels are reciprocal {i.e., hj^t = if 
the inter-user channel is not in outage, which happens with 
probability (1 — P e ), the four different packets transmitted to 
the BS are I\, I2, I\ ® I2 and 1\ Oh- Upon receiving two 
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Fig. 2. Two-user cooperative network with nonbiliary network coding, (a) 
Each user broadcasts its own information and (b) each user transmits a linear 
combination over GF(4) of all the available information packets. 



copies of the same message, the BS performs maximum ratio 
combining (MRC). In this case, an outage for Ii © I 2 occurs 

if (\hi,o. 2 \ 2 + \h 2 ,o, 2 \ 2 ) < 9 and has probability P x w \ 0. 

Without loss of generality, we analyze the outage probability 
of User 1 . The same result holds for User 2 due to symmetry. 
An outage event occurs when the direct transmission of T\ and 
one or two packets from I 2 and 1\ © I 2 cannot be decoded by 
the BS. In this case, the outage probability for User 1 is 

Pp,l «P e (P e +Pl). 

If the inter-user channel is in outage, which happens with 
probability P e , each user just retransmits its own information. 
Upon receiving two copies of the same information packet, 
the BS performs MRC, resulting in the following outage 
probability for User 1 J2] 

pi 

Pf 1 « — ■ 

/,i 2 

Considering all the possibilities for the inter-user channel, the 
outage probability of User 1 is then given by 

Po,l = (1 - Pe)Pp,l + PP/,1 « Pi 

If the inter-users channels are non-reciprocal, it can be 
shown that the outage probability of a 2-user binary network 
coded system is given by J4) 



Po 



BNC 



pi 



(6) 



which corresponds to a diversity order D = 2 according to 
(0. We can see that the diversity order obtained from (|6]l is 
not optimal, since the information of each user is transmitted 
through three independent paths and higher diversity order can 
be achieved, as will be explained later. By a similar analysis, 
the outage probability of the 2-user DAF scheme, presented 
in Fig. Q] considering non-reciprocal inter-user channels, can 
be shown to be 



Po.DAF ~ l-5-P e 



(7) 



C. Dynamic-Network Codes 

In @, (6), Xiao and Skoglund showed that the use of 
nonbinary network coding is necessary to achieve a higher 
diversity order, and then proposed the so-called DNC. A 
simple 2-user DNC scheme is presented in Fig. [2 where 
nonbinary coefficients are used. 

Assuming reciprocal inter-user channels and considering 
that they are not in outage, which happens with probability 



(1 - P e ), the BS receives I u I 2 , h ffl h and I x ffl 2I 2 . We 
can see that the BS is able to recover the original messages 
Ii and I 2 from any 2 out of the 4 received packets. For User 
1 (the same result holds for User 2 due to symmetry), an 
outage occurs when the direct packet I\ and at least 2 out of 
the 3 remaining packets cannot be decoded. This occurs with 
probability 

P p ,i=Pe (f)P?(l-P e ) + P2 «3P e 3 . 

However, when considering the link failure of inter-user 
channels, with probability P e a user cannot decode its partner's 
message. In this situation, it retransmits its own information. 
The BS performs MRC, and the outage probability for User 



1 is then P 



fA 
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13. Thus, the overall outage probability 



for User 1 is 

Po,l = PePfS + (1 - Pe)Pp.l « 3.5P e 3 . 

It is easy to see that the diversity order is D = 3. If the inter- 
user channels are non-reciprocal, then the outage probability 
was shown to be |5) 



P 



0,1 



4P? 



For M users, in the DNC scheme, each user transmits a 
fixed number of M — 1 nonbinary linear combinations in the 
cooperative phase. It is shown in [5| that the diversity order 
achieved by DNC is 



D DNC = 2M - 1, 



(8) 



however, with a low, fixed rate R = M/AI 2 = 1/M. 

The network code necessary to achieve the diversity order 
in ^ was proved to exist through an association between 
the DNC scheme and the single-source linear robust network 
codes iflOl Th.l 1]. The network code design proposed in J3 is 
to find the transfer matrix which is nonsingular for all possible 
error patterns in the inter-user channels. This network design 
method can have high complexity as M increases. 

In (6l, a simplified construction of DNCs was presented. 
This reduction in the design complexity is obtained at the 
expense of a poorer performance, which is negligible if the 
field size is as large as H 



q> 



M 2 - 1 
M-l 



(9) 



III. Motivation for this Work 



The multiple-source-one-destination network can be repre- 
sented by a transfer matrix which in turn can be seen as 
a generator matrix of a systematic linear block code. The 
generator matrix obtained from the 2-user system, illustrated 
in Fig. 12 is given by 



»DNC 



" 1 
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1 " 





1 


1 


2 



(10) 



In the DNC scheme (5), ]6], the diversity order is related 
to the minimum number of correctly received packets (or 
symbols) at the BS from which the information packets from 
all users can be recovered. A package which is not received 
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correctly may be thought of as an erasure, and is discarded by 
the receiver. The receiver's ability to recover the information 
packets from the non-erased packages is thus equivalent to 
the erasure correction capability of the associated (network) 
block code. It is well-known that the transmitted codeword of 
a linear block code with minimum Hamming distance D can 
be recovered if no more than D — 1 of its positions have been 
erased by the channel 11201 . The connection between these two 
problems establishes that the diversity order of the 2-user DNC 
system shown in Fig. [2] under the assumption of perfect inter- 
user channels, is equal to the minimum Hamming distance of 
the rate 2/4 block code with generator matrix given in dTUb . 

In general, for a rate k/n linear block code, the mini- 
mum Hamming distance is upper bounded by the Singleton 
bound EO): 

<4rin<n-fc + l. (11) 

As shown in lfl8l Chap. 11, Coroll.7], the Singleton bound can 
only be achieved if the alphabet size is large enough. For 
example, for a 4/8 block code, the Singleton bound gives 
dmin < 5. However, the maximum possible achieved minimum 
Hamming distance in GF(2) is 3. In GF(4), it is possible to 
achieve rf m ; n = 4. The upper bound rf m ; n = 5 is only possible 
to be achieved if the field size is at least 8 l2D . 

In the DNC scheme 0, 0, the overall rate is M/M 2 . 
From ( fTTT i, the diversity order is thus upper bounded by 
-Dmax.DNC = M 2 — M + 1. We observe that a system with 
rate R = aM/aM 2 (for a > 2) would have the same overall 
rate as the DNC scheme (R = 1/M), but the diversity upper 
bound would be increased to 



a{M 2 - M) + 1. 



(12) 



This motivates us to modify the DNC scheme accordingly 
and to propose an even more general scheme which is more 
flexible in terms of rate and diversity order. 

Due to inter-user channel outages, the upper bound pre- 
sented in (TT2l cannot be achieved. In Section [V] the discrep- 
ancy between the upper bound and the real diversity order 
obtained is quantified by the outage probability analysis. 

IV. Generalized Dynamic-Network Codes 

We begin the description of the proposed GDNC scheme 
by elaborating on the 2-user DNC scheme presented in Fig. [2] 
whose associated generator matrix is given in ( [Tol l. The new 
scheme is illustrated in Fig. [3] 

Each user broadcasts three packets (or symbols) of its 
own in the broadcast phase, and then each user transmits 
two nonbinary linear combinations (of the six previously 
broadcasted packets) over GF(q = 2 b ) in the cooperative 
phase, where b is an integer greater than zero. The receiver 
collects the 10 packets, which can be seen as a codeword of 
a systematic 6/10 linear block code. 

Without loss of generality, we analyze the outage probability 
for User 1 in the first TS. Let us assume for the moment 
that the inter-user channels are reciprocal. In this case, the 
probability that no errors occur in the inter-user channel is 
1 — P e . If User 2 can correctly decode -Zi(l), the message 
of User 1 in the time slot 1, then 7i(l) will be in outage at 
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(a) Broadcast phase 1. (b) Broadcast phase 2. (c) Broadcast phase 3. 




O' 

User 2 
(d) Cooperative phase 1. 
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User 2 
(e) Cooperative phase 2. 



Fig. 3. A rate 6/10 GDNC scheme with M = 2 users. The symbol 
EB m t i denotes the linear combination of all the available information packets 
performed by relay m in time slot t' of the cooperative phase. 



the BS if the direct packet and, in the worst case, the 4 parity 
packets containing ii(l) cannot be decoded by the BsQ, which 
happens with probability P p i w P e 5 . 

If User 2 cannot correctly decode h{l), which happens with 
probability P e , it will not be able to help User 1 by relaying 
-Zi(l). In this case, the BS receives only 3 packets containing 
2i(l) (the direct transmission plus 2 parity-checks sent by 
User 1 itself). Message ii(l) will then be in outage at the BS 
if the direct packet and, in the worst case, both of the parity- 
checks transmitted by User 1 cannot be decoded either. This 
happens with probability Pf i « P^. Therefore, considering 
reciprocal inter-user channels and all the outage patterns, the 
outage probability of message ii(l) is given by 



Po.l = PePf.l + (1 - Pe)P Pt l « Pe 



(13) 



We can show that the same result is obtained when the inter- 
user channels are non-reciprocal. 

We can see from ( fT3l that the diversity order achieved by 
the rate 6/10 GDNC scheme with M = 2 users presented in 
Fig. [3] is D = 4, which is higher than the one obtained by 
the rate 2/4 DNC scheme of 0, with M = 2 users in ©. 
Therefore, both the rate and the diversity order have increased. 

We can also verify that the outage probability is dominated 
by the term related to the inter-user channel being in outage, 
when User 2 cannot help User 1. Similarly, for a M-user 
network, we will see that when all the M — 1 inter-user 
channels are in outage is also the worst-case scenario for the 
outage probability of a given message. 

A. Multiple Users 

The generalization of the scheme presented in Fig. [3] to the 
case of M > 2 users is illustrated in Fig. [4] Ij (t) represents 

'if User 2 can correctly decode there are other situations yielding 

the same outage probability for message /i(l). But, for the moment being, 
we are not interested in the multiplicity of the outage patterns, but only in 
the diversity order. 
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Userl: /i(l) ••• Ji(fci) } Pi(l) ••• Pi(fe) 

User 2: J 2 (l) ■•• / 2 (*i) I P 2 (l) ••• ft(fe) 

UserM: J M (1) ■•• iaf(fci) [ Pm(1) ■•• P^(*a) 

time slots 

Fig. 4. GDNC scheme for a network with M users. 



the information (symbol or packet) transmitted by user j (j = 
1, • • • , M) in time slot t (t = 1, • • • , fci) of the broadcast phase 
(at left), and Pj (t') corresponds to the parity-check transmitted 
(symbol or packet) by user j in time slot t' (t' = 1, • • ■ ,k 2 ) 
of the cooperative phase (at right). 

Following Fig. |4] each user first broadcasts k\ independent 
information packets of its ow In the cooperative phase, each 
user transmits k 2 parity-check packets consisting of nonbinary 
linear combinations of its k\ own information packets and the 
kx(M-l) partners' information packets (if decoded correctly). 
If a user cannot correctly decode an information packet from 
one of its partners, this information packet is replaced by an 
all-zero packet in the formation of the linear combination. 
Herein we assume that, as in J5], the BS knows how each 
parity packet is formed. This may lead to some overhead, 
however, if I jit) is long enough, the overhead is negligible. 
Thus, the GDNC overall rate is given by: 

fciM ki 



Pgdnc — 



(14) 



k\M + k 2 M h + k 2 ' 

From (fTTT i and (TBI , we can see that the diversity order of the 
GDNC scheme is upper bounded by 



£>GDNC < k 2 M + 1. 



(15) 



However, due to the errors in inter-user channels, this upper 
bound cannot be achieved. In the next section, we find the 
achievable diversity order of the proposed scheme. 

V. Outage Probability and Diversity Order 

Denote by Dj t t Q {1, • • • , M} the index set corresponding 
to the users that correctly decoded I jit), the information 
packet of user j in time slot t in the broadcast phase. For 
convenience, include index j itself to Djj- The number of 
users in Dj :t is denoted by \Dj :t \- Let P e again be the 
outage probability of a single channel, and let Dj- t denote 
the complement set {1, . . . , M}\Dj t t, i.e., Dj :t contains the 
indices of the users which could not decode Ij(t) correctly. 

The probability of Djj is approximately p! 1 ^'*'. 

We define a new set T>j t t{I) as the set of all messages 
correctly decoded by the users in Dj.t in the broadcast phase, 
including Ij(t). There are at least \Dj :t iI)\ + \Dj,t\k 2 packets 
containing messages of Vj^il) (the \Dj^iI)\ messages them- 
selves from the systematic part and \Dj t \k 2 parities). Assum- 
ing that a well designed network code is used, the |X>j t (/)| 

2 These packets can be sent in any order. Recall that the channels are 
considered uncorrelated across both time and space. 



messages can be recovered from any \T>j_ t (I)\ received pack- 
ets^. For a fixed Dj t (which fixes Djj as well), an outage 
for Ij it) is declared only if the direct transmission Ij (t) and 
at least \Dj, t \k 2 out of the remaining |2?j, f (/)| + \Dj :t \k 2 — 1 
received packets are in outage, which occurs with probability 



|-Dj,t | &2 in outage 



P jiDj, t ) = P e 



V ht iI)\ + \D ht \k 2 -l 



\D jt t\k 2 



P! 



+ 



(l-Pe) 1 -!^.^)! ' 

V jtt (I)\ + \D jtt \k 2 - 1\ pl^WH^I^-i- 



V j>t (I)\ + \D j>t \k 2 -l 



(1 - Pe)° 



All the \V j:t (I)\ + \D jtt \k 2 -l in outage 

(16a) 



Pe 



\Vj, t (I)\ + \Dj, t \k 2 -l\ {Djt[k2 
\D h t\k2 



7(fci,*a, Dj, t )P t 



(M-\D jit \)k 2 + l 



(16b) 
(16c) 



where (?) is the binomial coefficient and 



Vj,t(I)\ + \Djj\k2 - 1 



\Dj,t\k 2 



is a positive integer representing the number (multiplicity) 
of outage patterns leading to that same probability. We can 
see that 7(fci, k 2 , Djj) increases as \Vj^iI)\ increases. In 
particular, since ki < [Dj t t{I)\ < Mk\, we must have that 



fc2 l<7(fci,fa,{l,...,M}\{j})<l ^ 

(17) 

It should be noted that ( TTSI l is the worst case, since the 
messages in Vj tt (I) (except Ij it)) might also be decoded by 
some users in Dj tt , reducing the outage probability of the 
messages in \T>jjiI)\. 

The outage probability is then given by 

Po-i = E Jf^'a - P^-^-^PoAD^t) (18a) 



Dj.t 



(M-\D ]tt \)k 2 + \D jtt \+l 



7(fci,fc 2 ,D Jit ) (18b) 



jik x ,k 2 , {1,..., M}\{j})P e M+k \ 



(18c) 
(18d) 



where P), iM (1 - p e )(M-i)-\Dj, t \ is the probability of \D j>t \ 
out of M — 1 inter-user channels in time slot t being in outage, 
and |-Dj,t|* corresponds to the -Dj,t value that results in the 
lowest exponent term in ( I18bb . which, for k 2 > 2, is |-Dj,t|* = 
M — 1. For k 2 — 1, the exponent in (II 8bb is equal to M + k 2 
independently of \Dj t \. We have proved the following result. 

3 For the moment being, we are just assuming the existence of such codes. 
We will prove their existence later. 
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Theorem 1. The diversity order of the GDNC scheme for an 
appropriately designed network code with sufficiently large 
field size is Dgdnc = M + k%- 

We can see that, when ki = 1 and k 2 = M — 1, the pro- 
posed scheme reduces to the DNC scheme Q, J6), with rate 
&i/(fci + k 2 ) = 1/M and diversity order M + k 2 = 2 A/ - 1, 
as mentioned at the end of Section IH-CI In particular, for 
ki = k 2 = 1 we have the 2-user DNC scheme shown in 
Fig. |2 The design approach in the current paper allows for 
different values of ki and k 2 , effecting different rates and 
different diversity orders. By varying ki and k 2 independently, 
we can obtain a tradeoff between rate and diversity. From ( TBI 
and Theorem [U we can see that an appropriate choice of k\ 
and k 2 can simultaneously improve rate and diversity order 
over the DNC scheme. 

VI. On the Network Code Design 

It was shown in Theorem Q] that the proposed GDNC 
scheme can achieve diversity order M + k 2 if the network 
code is appropriately designed. In this section, we design such 
network codes in the light of classical coding theory. First, we 
present some useful linear block codes properties. 

Let C be an (n, k, d m in) linear block code over GF(q) with 
a systematic generator matrix G given by 
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Recall from (fTTI) that the minimum distance of C is upper 
bounded by d m i n < n — k + 1. From now on, we assume that 
code C is MDS, i.e., that d m i n = n — k + 1. We will make 
use of the following results regarding the systematic generator 
matrix of MDS codes. 

Theorem 2. (OS Chap. 11, Th.8]) An (n, k, d min ) code C with 
generator matrix G = [I\P], as in d 1 9b . is MDS if and only if 
every square submatrix of P is nonsingular. 

The following lemma is a direct consequence of Theorem |2] 

Lemma 1. For any (n, k, d m i n ) MDS code C with generator 
matrix G = [I\P], the matrix P has no zero entry. 

Proof: From Theorem any one-by-one submatrix of P 
is nonsingular for MDS codes. Alternatively, any row of G is 
a codeword, and must have Hamming weight at least n—k+1. 

□ 

A. A Theory of Faulty Generator Matrices 

In this section, we study some properties of new block 
codes obtained from systematic MDS codes by zeroing some 
entries of its generator matrix, and we refer to the obtained 
generator matrix as a faulty generator matrix. In particular, we 
are interested in the minimum Hamming distance or the least 



erasure-correcting capability of MDS codes that can be assured 
when these codes are encoded by a faulty generator matrix. We 
will see that the theory of faulty generator matrices we develop 
next is related to and will shed light on our main problem: the 
design of network codes to maximize the diversity order of 
cooperative wireless communications assuming possible link 
failures. 

When a set of entries of the generator matrix P of the MDS 
code C is replaced by zeros, the new code C produced is no 
longer MDS. Let A = {(ai, &i), (a 2 , b 2 ), ...} be a subset of 
entries of P that are replaced by zeros, where (a, b) is the 
entry in row a and column b of matrix P. We call A a fault. 
Let A = {A , Ai, . . . , be a collection of / faults. We 

consider that two different faults Ai and Aj cannot contain 
a common entry of the matrix P. That is, Ai n Aj = 0. It 
is also considered that every fault has fixed cardinality, i.e., 
\Ai\ = \A\ Vi. Although these considerations may look too 
restrictive, we will see that they are sufficient for the purposes 
of solving our main problem. 

Let x = (XOi • • • j Xf-i) t> e tne binary indicator vector 
associated with the occurrence of faults, where 

Xi = 1 if Ai occurs 
Xi =0 if Ai does not occur. 

For a nonnegative integer i, let b(i) denote the binary (vector) 
representation of i. We denote the collection of all possible 
combinations of faults by B = {B^q), -B&(i), B^M-i)}, 
with B x = {(a,b)\(a,b) G U» :Xi =i Each event B X' 
which consists of the occurrence of |i? x |/|A| faults, where 
| ■ | stands for the cardinality of the set, gives rise to a 
new generator matrix of a block code C(B X ) with minimum 
Hamming distance d m i n (C(B x )). We also define the minimum 
composite distance of the code C{B X ) as 

<C?(C(B X )) = d mm (C(B x )) + \B X \/\A\, (20) 

which is composed of its minimum Hamming distance 
d m in{C(B x )) plus a "compensation" term related to the num- 
ber of faults in the combination B x . 

A parameter of fundamental importance to indicate the 
performance of a MDS code subject to a set of faulty generator 
matrices is the least minimum composite distance of any 
possible combination of faults. In particular, given a MDS 
code C and a collection of / faults A, this distance is defined 
as: 

QC, A) 4 min {dZl (C(B X ))} , (21) 

B X £B 

where, when there is no confusion, it is simply called minimum 
composite distance. 

We next prove a sequence of results which constitute a 
theory of faulty generator matrices. 

Lemma 2. The weight between two codewords y and z is 
lower bounded by the difference between their individual 
weights, i.e., w(y EH z) > w(y) — w(z). 

Proof: From the triangular inequality, w(y' EH z') < 
w(y') + w(z'). If we set y' = ySz and z' = Bz, where Bz is 
the additive inverse of z, we obtain w(y EH z) > w(y) — w(z), 
which proves the lemma. □ 
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Theorem 3. Let C be an (n, k, d m i n ) MDS code with system- 
atic generator matrix G = [I\P]. The replacement by zeros of 
5 entries in any single row of the matrix P gives rise to an 
(n, k, d m i n — 6) code C. 

The proof is provided in Appendix lAl 

Lemma 3. Let C be an (n, k,n — k — S + 1) code (with 
generator matrix G' = [I\P']) obtained from an (n, k, n—k+1) 
MDS code C according to Theorem \3\ The replacement by 
zeros of additional entries in any of the 8 columns ofP' where 
zeros have previously been inserted gives rise to a code C" 
with the same minimum Hamming distance as C 

Proof: Each codeword of C" can be obtained from a 
codeword of C by changing 6 positions, which are fixed for 
all codewords. Therefore, we must have that d m i n (C") > 
d m in(C) — S. Since g' it (see Appendix lAb. a codeword of C 
of Hamming weight n — k — 6 + 1, is also a codeword of C", 
the Lemma is proved. □ 

Lemma 4. If C is an (n, k, d m i n ) MDS code with systematic 
generator matrix G = [I\P], the replacement by zeros of 5 
arbitrary entries of its matrix P generates a code C with 
minimum Hamming distance d m i n (C') > d m i n — S. 

Proof: Similar to Lemma [3] with equality if the new 
code C has a codeword with Hamming weight d m i n — 6. For 
example, if all the S replacements occur in a single row of P. 

□ 

According to Lemma @] a single fault A (confined to the 
matrix P) can reduce the minimum Hamming distance of 
the original MDS code by at most \A\. On the other hand, 
Theorem [3] tells us that the worst scenario for a single fault 
A is when all of its entries belong to a single row of P. The 
last result of this section refers to the worst scenario for a 
collection of faults. 

Lemma 5. Let C be an (n, k, d m i n ) MDS code with systematic 
generator matrix G = [1\P]. Consider a collection A = 
{Aq, A\ 1 . . . , Af-i) of f faults, where \Ai\ = \A\ > 1, Vi. 
If all the faults occur in a single row of P, then the minimum 
composite distance of C subject to A is a strictly decreasing 
function of the number of faults. 

Proof: From Theorem [3] each such a fault reduces 
the minimum Hamming distance by exactly \A\ > 1, and 
contributes with only 1 to the "compensation" term of the 
minimum composite distance. In the special case where 
| A | = 1, the minimum composite distance of C is constant 
independently of the number of faults. □ 

B. Maximum-Diversity GDNC 

Recall that our main problem is to design a linear network 
code for the GDNC scheme that, given all the inter-user 
error patterns and taking into account their probabilities of 
occurrence, assures that the diversity order is M + k%. The 
solution to this problem is our main result, presented next. 

Theorem 4. If a systematic generator matrix of a MDS code 
C with minimum Hamming distance d m j„ = Mk2 + 1 is used 



as a transfer matrix of the GDNC scheme, then the diversity 
order D CDNC = M + &2 is guaranteed. 

The proof of this theorem is provided in Appendix [B] 

Theorem 5. If a generator matrix of a non-MDS code is used 
as the transfer matrix in the GDNC scheme, the maximum 
diversity order presented in Theorem\l\is not guaranteed. 

The proof is provided in Appendix ICl 

Now that we have proved that MDS codes are sufficient to 
achieve the maximum diversity order of the GDNC scheme, 
the network code design problem is reduced to a MDS code 
design, without the hassle of testing a large number of network 
codes, taking into account all the inter-user channel error 
patterns, and their impact on the system's diversity. The well- 
known class of MDS codes called Reed-Solomon (RS) codes 
l20l Chap. 3] lfl8l Chap. 10], for example, may be used. In the 
destination node, the users' message packets can be recovered 
with some classical decoder for RS codes, e.g. l22l . l23l . 

The generator matrix of an (n, k,n — k + 1) RS (or 
extended RS) code over GF(g) in systematic form can be 
easily obtained with the software application SAGE l24l . 
entering the following code: 

F . <alpha> = GF (q) 

C = ReedSolomonCode (n, k, F) 

G = C . gen_mat ( ) 

show (G . echelon_f orm ( ) ) 

It should be noted that the field size q must be sufficiently 
large lfl8l Chap. 11, Coroll.7]. In this regard, we have the 
following result. 

Theorem 6. There exists an MDS code (n = M(kx+k2),k = 
Mki,d m i n = Mk2 + 1) over GF(q) (where q is a power of 
prune) that can serve a GDNC system achieving maximum 
diversity order M + ki if q > M{k\ + £2). 

The proof is provided in the Appendix [D] 

We can now compare the field size necessary for the full 
diversity order to be achieved in the GDNC scheme (according 
to the bound presented in Theorem [6]l and the simplified DNC 
scheme (as proved in J6) and shown in (0). This comparison 
is exemplified in Table U 

TABLE I 

Field size necessary for the diversity order achieved by the 

SCHEMES GDNC (WITH ki = 1 E k 2 = M — 1) AND DNC TO BE 
GUARANTEED. 



M 


2 


3 


4 


5 


</ min (GDNC) 


4 


9 


16 


25 


5min (DNC) 


3 


28 


455 


10626 



We emphasize that the values shown in Table U regarding 
the DNC scheme correspond to the finite field size values 
necessary for the diversity order 2M — 1 to be achieved, 
according to ©. However, since the GDNC scheme with 
parameters ki = 1 and k% = M — 1 reduces to the DNC 
scheme, we can see that our result in Theorem [6] represents a 
much tighter lower bound on the field size. 
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Tables [TT] and [TTT] list the parity matrix P described in (fT9] l 
of some MDS codes obtained from RS codes (as explained in 
the proof of Theorem [6} which can be used as the GDNC's 
network transfer matrix, for M = 2 and M = 3 users, 
respectively. These matrices were obtained using SAGE 



link level with which it is possible to recover the transmitted 



TABLE II 

Network codes obtained from RS codes for M ■■ 



: 2 USERS. 



kl 


k 2 


Code rate 


Field size q 


Parity matrix P 


1 


1 


2/4 


4 


"3 2" 
.2 3. 


1 


2 


2/6 


8 




"3 5 17" 
.2 4 3 6. 




2 


1 


4/6 


8 


Obtained from the 4/8 code by 
puncturing the last 2 columns 


2 


2 


4/8 


8 




"3 7 3 6" 

5 7 7 4 
2 4 6 1 
.5 5 3 2. 




2 


3 


4/10 


16 


Obtained from the 4/12 code by 
puncturing the last 2 columns 


2 


4 


4/12 


16 




"14 5 7 10 4 10 8 9 " 

12 4 14 5 12 12 14 7 

13 6 12 15 8 1 4 10 
.14 64 1 1 6 3 5_ 




3 


2 


6/10 


16 


Obtained from the 6/12 code by 
puncturing the last 2 columns 


3 


3 


6/12 


16 




-11 2 4 6 14 12" 

1 11 13 10 14 10 

2 4 2 10 5 9 
6 13 12 11 8 12 
4 12 12 2 6 6 

-11 13 10 14 10 4 - 





TABLE III 

Network codes obtained from RS codes for M ■ 



■ 3 USERS. 



fcl 


fc 2 


Code rate 


Field size q 


Parity matrix P 


1 


1 


3/6 


8 




"4 1 5' 

3 13 
.6 1 7. 




2 


1 


6/9 


16 


Obtained from the 6/12 code by 
puncturing the last 3 columns 


1 


2 


3/9 


16 




"15 11 1 14 5 11" 

8 5 1 8 13 4 
. 6 15 1 7 9 14. 




2 


2 


6/12 


16 




-11 2 4 6 14 12" 

1 11 13 10 14 10 

2 4 2 10 5 9 
6 13 12 11 8 12 
4 12 12 2 6 6 

-11 13 10 14 10 4 - 





VII. Simulation Results 

In order to support the outage probability analysis of the 
previous section, we have performed some computer simula- 
tions. The frame error rate (FER) was simulated and plotted 
against the SNR. In our simulations, we adopt r = 0.5 bits per 
channel usfl We assume that there exists a channel code at the 

4 It should be noted that the choice of the information rate r is irrelevant 
for the purpose of obtaining the diversity order. 



packet if \hj^ t t 



> 



if I^",mI 2 



< g, an outage in the link 



j — > i is declared, and the generator matrix G will have a set 
of its elements replaced by zeros. The number of erroneous 
frames is computed based on the rank of the resulting faulty 
generator matrix. The network codes used in the simulations 
were obtained with the method presented in Section [VT] The 
analytical outage probabilities obtained in this paper were also 
plotted (dashed lines). 

Fig.[5]shows the FER versus SNR for a 2-user network, con- 
sidering the BNC scheme JU, the DNC scheme 0, @] (Fig. 0) 
over GF(4), and the proposed GDNC with fci = k 2 = 2 
over GF(8), with the generator matrix obtained from Table [TT] 
All of them have the same rate equal to 1/2. Regarding the 
analytical GDNC outage probability, we consider an upper 
bound which was obtained by substituting the upper bound of 
7(fcij ^2, Dj_ t ) in (l% into d 1 8db - As expected, the proposed 
scheme achieves a higher diversity order than the two other 
schemes with the same rate. The SNR gap between the 
analytical and simulated curves are due to the Gaussian input 
assumption made in (f2j) and the use of the aforementioned 
upper bound. However, we can see that the simulated diversity 
order (curve slope) matches the analytical one. 



w 10 




SNR (dB) 

Fig. 5. FER versus SNR (dB) for a 2-user system with rate R = 1/2, 
considering the BNC scheme, DNC scheme (in GF(4)) and the proposed 
GDNC scheme (with fci = &2 = 2 and in GF(8), according to Table Hit. 

It should be mentioned that, since the maximum diversity 
order 3 of the DNC scheme 151, l6l is already achieved 
with the field GF(4), increasing the field size would not 
bring any advantages in this case. On the other hand, as 
discussed in Section [III] for the proposed GDNC used in 
our simulations, which corresponds to a rate 4/8 block code, 
the field size 8 was necessary for d m i n = 5. It should also 
be pointed out that, in the proposed scheme, only the linear 
combinations are performed over GF(q) prior to modulation 
and transmission; the modulation scheme and the transmission 
mode are irrelevant to the functioning of GDNC, and can be 
the same for all schemes. 

Fig. [6] shows the FER performance for the proposed GDNC 
scheme illustrated in Fig. [5] as well as for two other GDNC 
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schemes with non-MDS codes (see Table |IV| ) with the same 
parameters k% = k2 = 2 and rate 1/2, but with respective 
block codes of minimum Hamming distances 3 and 4, whose 
generator matrices have been randomly generated. 



DNC (simulation) 
o- ■ DNC (outage prob) 
^ — GDNC, d min =3 (simulation) 

GDNC, d min =4 (simulation) 
* — GDNC, d m|n =5 (simulation) 
*- GDNC (outage prob) 




Fig. 6. FER versus SNR (dB) for a 2-user system with rate R = 1/2, 
considering the DNC scheme (in GF(4)) and the proposed GDNC scheme 
(with fci = &2 = 2 and in GF(8), according to Table ITT1 and HVt . 

We can see that, with d m i n = 3, the diversity order is only 
3, while with d m i n = 4 the diversity order is close to the one 
obtained with the proposed code (d m in = 5). In this case, a 
quasi-MDS code achieves the maximum diversity. 

TABLE IV 

Network codes obtained from codes with minimum distance 3 

AND 4. 



Fig. 


^min 


Parity matrix P 


[6] and E] 


3 




"0 5 5 3" 
4 7 7 1 
7 7 7 1 

.7 2 3 7. 




□ 


4 




"4 2 5 3" 
5 4 6 4 
7 6 14 

.113 1. 




El 


4 




"0 5 5 5" 
7 4 6 3 
7 4 7 

.2 7 6 0. 





However, in Fig. [7] a different code, also with d m i n = 4 
(given in Table HVl i. is shown to fail to achieve the maximum 
diversity order. Thus, we can see that by using a generator 
matrix of a quasi-MDS block code as the network transfer 
matrix of the GDNC scheme, the diversity order M + k^ is 
not guaranteed. This is in accordance with Theorem [5] 

VIII. Conclusions and Final Comments 

In this work, we have proposed a generalization of the 
network coding method in (5). The aim of the generalized 
dynamic-network code (GDNC) is to increase the diversity 
order of cooperative wireless communications systems without 



-0— DNC (simulation) 
-o- DNC (outage prob) 
-0 — GDNC, d min =3 (simulation) 



. GDNC. d =5 (simulation) 

mm * ' 

GDNC (outage prob) 




Fig. 7. FER versus SNR (dB) for a 2-user system with rate R = 1/2, 
considering DNC scheme (in GF(4)) and the proposed GDNC scheme (with 
fci = fc2 = 2 and in GF(8), according to Table ITll and UVt . 



sacrificing the system's rate, what appears to be a drawback 
of the original method. We have shown that the problem of 
designing network codes that maximize the diversity order is 
related to that of designing optimal linear block codes, in the 
Hamming sense, over a nonbinary finite field. 

An outage probability analysis was presented and computer 
simulations supported the analytical results. The two design 
parameters fci and &2 of GDNC may be varied to produce a 
wide range of rates and diversity orders, and offer a tradeoff 
between rate and diversity order (when the value of &2 is 
changed), or even between rate and decoding latency (when 
ki is changed), as compared to the baseline system. 

Of course, since the GDNC scheme improves the rate- 
diversity trade-off by coding across more time slots, the 
decoding complexity and the delay are increased. However, 
for relatively small number of users and small values of fci 
and fc2, the decoding complexity and the delay are manageable 
as efficient decoding methods for RS codes do exist l22l . l23l . 
On the other hand, for a large number of users and/or when 
the parameters k\ and k^ are very large, the dimensions of the 
generator matrix of the block code associated with the GDNC 
scheme may be too large for a practical implementation of the 
encoders and decoders. Moreover, when the diversity order is 
already high, a further increase does not bring much additional 
benefit. In this case, the generator matrix can be made sparse 
and the overall system can be seen as a generalization of 
the adaptive network coded cooperation (ANCC) method 
proposed in l25l . The connection between GDNC and ANCC 
is currently being investigated. 

Appendix A 
Proof of TheoremO 

Proof: Let G = { gl , . . . ,gj and G' = {g[, . . . ,g' k } be 
the generator matrices of C and C , respectively, with g, and 
g' ; being the ith row of G and G', respectively. Without loss of 
generality, we assume that 5 elements of the parity part of gj , 
are replaced by zeros for some i* G {!,...,&}. For a nonzero 
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information message u = {ui, . . . ,Uk}, the Hamming weight 
of the corresponding nonzero codeword v' is 

w(V) = w(uG') 

= u;(«igi ffl. . .ffl tii.gj. ffl. . .ffl u k g' k ) 

( => w{u lgl ffl. . .ffl «i.gj. ffl. . .ffl Ufcg fc ffl« i .g i . Bui-gj.) 
5, 



the one that generates the least composite minimum distance, 
which is then given by 



(b) 

> w{Ul%i ffl. 

> (n - k + 1) 



where (a) follows from the fact that = g' ; Vi 7^ i*, (b) 
follows from Lemma [2] and (c) follows from the Singleton 
lower bound. Without loss of generality, we assume that 5 
elements of the parity part of <?.;» are replaced by zeros, giving 
rise to <?,■», a codeword of C of weight n — k + 1 — 5. The 
theorem is proved. □ 

Appendix B 
Proof of Theorem|4] 

Proof: Relating our problem with the one in Section 
IVI-AI one fault corresponds to one inter-user channel being in 
outage. When that happens, the user's receiver cannot correctly 
decode its partner's information, so, when forming k 2 linear 
combinations to generate its k 2 parity-check packets, this 
user replaces this erroneous packet with an all-zero packet, 
or, equivalently, sets to zero the k 2 coefficients associated 
with this partner. This amounts to replacing by zeros the k 2 
corresponding entries (in same row) of the parity matrix P, 
i.e., \A\ = k 2 . Since each user knows its own information, k 2 
entries in each row of P are immune to faults, while the other 
k 2 (M — 1) entries are subject to faults. In the worst scenario, 
when all the possible faults happen, the generator matrix takes 
the form 



Pi 



Pi 



M 



(23) 



where the (k\ x k 2 ) submatrix Pj contains the immune entries 
associated with user i. From Theorem [2] we know that every 
submatrix of Pi is nonsingular. Thus, the least minimum 
Hamming distance of a block code obtained from the original 
MDS code C due to the occurrence of faults is k% + 1. 
Nevertheless, the same minimum Hamming distance can be 
achieved with a much lower number of faults. From Lemma[3] 
and Theorem [3] we can see that the minimum number of faults 
for a code with minimum distance k 2 + 1 is M — 1, when all 
of these faults occur in the same row of P, for example. 

For all the possible minimum distances in the range k 2 + 
1 < dmin < AIk 2 + 1, a sufficient condition for the worst 
possible scenario (the lowest number of faults that result in 
this minimum distance) is when all the faults are located in the 
same row of P, according to Lemma [4] Thus, from Lemma [5] 
we can see that the larger the number of faults in a given 
row (and consequently the lower the minimum distance of the 
resulting code), the lower the composite minimum distance. 
This assures that the code with minimum distance k 2 + 1 is 



jcomp 



min {d 

B X €B 



(C(B X )) + \B X \/\A\} 

(k 2 + 1) + (M - 1) 
M + k 2 



(24) 



It is easy to see the connection between the two terms in the 
composed minimum distance and the exponents of P e in d!6cb 
and d 1 8ab . The proof is now complete. □ 



Appendix C 
Proof of Theorem[5] 

Proof: Let C be an (n, k, d m i n ) MDS code with generator 
matrix G = [I|P] . Replacing any set of entries of its parity 
matrix P by zeros gives rise to an (n, k, d' min ) non-MDS code 
C with d' min < d m i n , according to Theorem[3] Assume that C 
is used as the non-MDS network code for the GDNC scheme. 
If any of the entries replaced by zero is one of the entries that 
are immune to fault, we can easily show that there is a pattern 
of M— 1 faults that generates a codeword with weight < k 2 + l. 
Thus, the composite minimum distance is < M + k 2 , 

according to d2TT i. □ 

Appendix D 
Proof of Theorem[6] 

Proof: We want to find an (n = M(ki + k 2 ),k = 
Mki,d m in = Mk 2 + 1) MDS code over GF(q). Let us 
start by choosing a field size q > M(ki + k 2 ) as a power 
of a prime number. From JT8), there exists an (n' = q — 
l,k' = Mk x ,d' mm = q - Mkx) RS code over GF(g). 
If q = M(ki + k 2 ) + 1, simply use the RS code as the 
MDS code. If q = M(ki + k 2 ), we can add a parity-check 
symbol and have an extended (q,Mk\,q — Mk\ + 1) RS 
code over GF(g), which can be used as the MDS code. If 
q > M(k\ + k 2 ) + 1, then we can puncture the RS code 
in exactly q — M{k\ + k 2 ) — 1 positions, resulting in an 
(n = M(ki + k 2 ), k = Mki, d m i n = Mk 2 + 1) code, which 
is MDS. □ 
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